Coherent transport in disordered metals: zero dimensional limit 
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We consider non-equilibrium transport in disordered conductors. We calculate the interaction 
correction to the current for a short wire connected to electron reservoirs by resistive interfaces. 
In the absence of charging effects we find a universal current-voltage-characteristics. The relevance 
of our calculation for existing experiments is discussed as well as the connection with alternative 
theoretical approaches. 
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In recent years considerable attention has been de- 
voted to the effects of the Coulomb interaction on the 
transport properties of small structures, like thin diffu- 
sive films and wires [1, ^, ^, tunnel junctions |H, ^ 0, 
and quantum dots One interesting issue concerns 
the way an applied bias voltage affects the interaction 
corrections to the electrical conductivity. In diffusive 
metals these corrections arise from the combination of 
the electron-electron and impurity scattering and yield 
well known singularities at low temperature It has 
been shown that a finite voltage or, more in general, a 
non-equilibrium situation leads to a suppression of these 
singularities 0, |ll|, |l|, . In particular, in |lo[ |l| 
non-equilibrium transport in a short wire connected to 
electrical reservoirs by ideal interfaces has been consid- 
ered. However, in actual experiments the interfaces need 
not be ideal. Recently Weber et al. |15 investigated 
experimentally the non-equilibrium transport through a 
metallic nano-scale bridge. Both in ||l^ and in it has 
been suggested that the Coulomb interaction effects are 
responsible for the observed temperature dependence of 
the conductance and the current-voltage-characteristics. 
Whereas found an agreement between theory and 
experiment starting from a tunneling Hamiltonian, fl^ 
pointed out that the experimental data agree with what 
they expect for a diffusive conductor. In this paper we 
develop a formalism in which both the resistive behavior 
due to the interfaces and due to the diffusive wire region 
are treated on the same footing. From our results we 
conclude that the main resistive behavior in jist occurs 
at the interfaces. 

To begin with we recall the classical description of elec- 
trical transport through structures consisting of both in- 
terface barriers and diffusive regions. To be definite we 
consider a system made by a diffusive wire of length L 
which is attached to the reservoirs by two interface barri- 
ers. We study the system in a non-equilibrium situation 
with an applied voltage Vi — Vr — V where the sub- 
scripts I and r indicate the left and right reservoirs, re- 
spectively. The classical resistance of the structure is the 
sum of the wire resistance and the interface resistances 
-Rtot = -Rwiro + Ri+ Rr, so that the current as a function 
of voltage is / = V/Rtot- The microscopic calculations 



are conveniently carried out by using the Keldysh formal- 
ism |]l7| . In a disordered system the Keldysh component 
of the Green function reads as 



Gf (x,x) 



F,(x)[Gf(x,x) 
-27riiVoFe(x). 



Gf(x,x)] 



(1) 
(2) 



The first line is an exact relation and defines the distribu- 
tion function F. In the second line it is assumed that the 
density of states is a position independent constant. The 
current flowing in the wire or through the boundaries is 
then given by 

eDNoA Jded^F,{x) (3) 

// = er,^iVo J de[F,(0)-F,'] (4) 

Ir = eTrANo j de[F: - F,{L)l (5) 

where D is the diffusion constant, A the cross section, 
Ti^r are the interface transparencies, x — ... L is the 
position along the wire. The reservoirs are assumed to 
be in thermal equilibrium, with the distribution function 
given by F^^^ = tanh [(e — eVi_r)/2T] . The boundary con- 
ditions 



Dd^F\,^o 
Dd.FL^L 



~ri[F'~F{0)] 

-r,[F(i)-Fn 



(6) 
(7) 



guarantee current conservation at the interfaces. Fur- 
thermore we assume that the wire is so short that we can 
neglect inelastic scattering. In this case the kinetic equa- 
tion for the distribution function inside the wire becomes 
—Dd^F — 0. After solving this equation with the appro- 
priate boundary conditions one finally finds the current 
as / = V/{Ri + i?wiro + Rr) with 



Rr 



Ri = 



l/{2e^ANoTi) 
L/{2e^ANoD) 
l/{2e^ANQTr) 



(8) 
(9) 
(10) 



as one expects for three resistors in series. 
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We now discuss the quantum correction to the current. 
FoUowing we divide the quantum correction to the 
current in two contributions: The first contribution has 
the meaning of a correction to the conductance, SI^^"^ — 
V6G. The second one is due to the quantum correction 
of the distribution function and can be interpreted as the 
redistribution of the voltage along the interfaces and the 
wire, SI^^"^ = G5V. Both terms are necessary in order 
to ensure current conservation. By exploiting current 
conservation in the structure, 6Ii = (5 /wire — Sir, and 
fixing the voltage drop over the whole system to V, so 
that the sum SVi + 6Vr + (5Kviro is zero, it is possible to 
eliminate J/'^^^ from the above equations to get for the 
correction to the current : 



SI = 



RiSI, 



(1) 



(1) 



(1) 



Rl + Rr + i?„ 



(11) 



To proceed further we need the explicit form of SIi^J^ 
For an interface attached to an ideal lead on one side the 
quantum correction to the current is controlled by the 
correction to the density of states on the other side. 
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(1) 



eATi / deSN{e,0)[F,{0)-Fi] 



(12) 



64^^ = eATr / deSN{e,L)[Fl' - F,{L)]. (13) 



The limit where both sides of the interface are in thermal 
equilibrium has been studied many times in the literature 
[^j, . Out of equilibrium we obtain the density of states 
correction as 



6N{e,x) = -No I ^S{x,x) 



(14) 

S{x,x) = Im J dxi 

xF^^^{xi)p^{x,xi)^^{xi,x), (15) 

where pi^{x,xi) describes the spreading of a charge in- 
jected into the system at xi; it satisfies the equation 



(— icj — Dd^)puj{x, x') — eS{x — x'). 



(16) 



The quantity ^i^{xi,x) is the electrical potential at xi of 
a charge that has been injected at a;. It is given by the 
product of the dynamically screened Coulomb interaction 
with the diffusion propagator 

e^$i^(a;,a;') = dxiVuj{x,xi)p^{xi,x'). (17) 



p and $ depend on the details of the device under con- 
sideration and we will come back to them below. 

The expression for the correction to the current in the 
wire has been obtained diagrammatically in |l4[ and is 
given by 



«lirc(^) - ^eDN.A I de 



+2eDNoA / de 



^d,[F,{x)Six,x)] (18) 

'^F^{x)da:^S{x,Xi)\ooi=x- 



In general SI^^j.^{x) depends on the position x. Eq. ( pT| ) 
however is constructed in such a way that the spacial 



average J/ ^,"^^^ — L ^ J dxi5/^^;^g (x) has to be inserted. 



Eqs.( |12| , |l3yi^ ) allow to calculate the quantum correc- 
tions far from thermal equilibrium for an arbitrary ge- 
ometry. One observes that the only ingredients are the 
position dependent distribution function F, the charge 
density p, and the field We have already discussed the 
distribution function F including the relevant boundary 
conditions in Eqs. (|^), (^ and below. Inside a diffusive 
wire the charge density p satisfies Eq. ([1^). At the bound- 
aries with the left and right reservoir one may derive the 
matching conditions 



DdxPu,{x,x')\ x=0 

DdxPuj{x,x')\x=L 



-TrPuiL^x'), 



(19) 



to be compared with Eqs. (^) and (0). A careful analysis 
is also required for the field ^^^{x,x'). In the case of good 
metallic screening an injected charge is almost instantly 
screened so that the wire will be electrically neutral with 
the exception of a thin surface layer. In this case one has 
inside the wire 



ad'^^^{x, x') = e5{x — x'), 



(20) 



where a = 2e^DNQA is the conductivity. In the absence 
of surface charges the boundary conditions for $ and p 
are identical. In the presence of these charges, however, 
this is not the case. 

Many special cases where the general formalism dis- 
cussed above applies have been discussed in the litera- 
ture. For example standard Coulomb blockade physics 
is found when the resistance of the system is dominated 
by one of the interfaces |l^. In this limit a non- 
equilibrium analysis is not necessary since the distribu- 
tion function has the equilibrium form on both sides of 
the interface. The limit of two highly resistive interfaces 
has been studied in ||l5|, 20 and the opposite limit of a 



resistive wire and no interface barriers has been discussed 
in ||l^, Q . In the following we concentrate on a short re- 
sistive wire with interfaces assuming temperatures of the 
order of and lower than the Thouless energy, hD/L^. 

We start by expanding the charge density p in diffusive 
modes, 



Puj{x,x') = 



fn{x)U{x') 



(21) 



where the (normalized) functions fn{x) are obtained by 
the eigenvalue equation 



Ddlfn{x) = 7„/„(x). 



(22) 



In the zero dimensional limit we approximate the sum 
in Eq. ( ^l] ) by retaining only the eigenmode with the 
lowest energy, i.e. Pu,{x,x') efo{x)fo{x')/{~iuj + 70). 
This approximation is justified when the energy scales 
related to the temperature and to the voltage remain 
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FIG. 1: Temperature dependence of the conductance for dif- 
ferent values of the bias voltage. dl/dV is in units of /h 
and must be multiplied with the non-universal number A de- 
fined in the Eq. (feq). 7o is the energy of the lowest diffusive 
mode; in the absence of interface barriers 70 — n'^hD/L^, for 
strongly resistive interfaces 70 0. 



below the energy of the second lowest diffusive mode. Let 
us for the moment ignore charging effects. Then the field 
^ui{x, x') is frequency independent and one observes that 
the frequency dependent factors in all the contributions 
to the current are identical, (S/(T, V) ~ F^_^F^' /{—iu} + 
7o). The explicit result reads 



SI = -A— 
2n 



pOO 


ttT 


/ d?7e-'^'"' 
'0 


sinh(7rT?7) 



sin(ey?7)(23) 



where only the dimensionless number A and the quan- 
tity 7o depend on the details of the system under con- 
sideration. Notice that the integral has to be cut off 
at short times in order to avoid a logarithmic diver- 
gence. Let us first discuss the temperature and volt- 
age dependence of SI, before determining A and 70 ex- 
plicitly in the two limits of perfectly transparent inter- 
faces and for interfaces with low transparency. Figure 
|l| shows dI/dV/{Ae'^/2TT) as a function of temperature, 
the classical conductance has been subtracted. At high 
temperature there is a logarithmic behavior, which sat- 
urates below Tsat ~ max(7o,ey). Figure shows the 
voltage dependence of the conductance. Note that the 
linear conductance has been subtracted. For 70 <C T the 
conductance scales with voltage over temperature, while 
when 7o is large the relevant scale for conductance vari- 
ations is 7o. 

How large are the amplitude A and the energy of the 
lowest diffusion mode 70? In the case of two well trans- 
mitting interfaces, Gwirc ^ Gi,Gr, the eigenfunction of 
Eq. ( p^ ) with the lowest eigenvalue is 



/o(x) ^dj- sin(7ra;/L), 70 = n^D/L^ (24) 

The distribution function and the potential $ are deter- 
mined as 



FIG. 2: Voltage dependence of the conductance. In all curves 
the linear conductance has been subtracted. 



^u^ix,x') = 



[L ~ x')x/L'^ X < x' 
{L - x)x' /L"^ X > x' 



(26) 



and the amplitude of the correction to the current is 
found to be A = 64/7r'' - A/tt"^ « 0.25. In the oppo- 
site limit, Gwirc ^ Gi,Gr the eigenvalue equation ( [2^ ) 
may be solved perturbatively in the barrier transparency 
and one obtains 



fo{x) = l/VL,7o = (r, + r,)/L (27) 

Fix) = (r,^^, + r,F,)/(r, + r,) (28) 

$^(a:,x') = e/{Gi + Gr), (29) 

which leads to ^ = 2TiTr/iTi + Tr)'^. 

It is useful at this point to briefly discuss how the 
charging effects modify the above results. In the case 
of highly transmitting interfaces, we assume that the ac- 
cumulated charge is proportional to the field as p{x) = 
G<I>(a;), where G is a capacitance per unit length. This 
leads to a diffusion equation for the field 



— iwG$cj(x, x') — ad^^uiix, x) = eS{x — x) 
with the solution 

sin(n7ra;/L) sin(n7ra;'/L) 



$^(a:,a:') -2e^. 



-ioj{CL) -f {-KnfG^ 



(30) 



(31) 



F{x) = [{L- x)Fi+xFr]/L 



(25) 



and the correction to the current is modified according 
to 

SI = -— dyye-"^" ■ u ^ ^^^i^^^) 
Ztt Jq [smh(7rJ r/) J 

X - exp[-r;(7rn)Vi?(GL)]}. (32) 

n 

The numbers An depend on the wave-function of the dif- 
fusive mode n and R — G^^^^. Notice that although the 
charge density pui{x,x') is zero dimensional many diffu- 
sive modes have to be taken into account in the field 
(j)ij{x,x'). For a small capacitance G the charging sim- 
ply cures the short time divergence in the integral in 
Eq. (|2^). For a larger capacitance the /-^-characteristics 
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is no longer a universal function. For the system with 
poorly transmitting interfaces, the field in Eq. ( p9|) has to 
be replaced by ^uj{x,x') = e/{—iujC + Gi + Gr), where 
we denoted the capacitance of the system by C. The 



modification to the current is analogous to Eq. (|32|). 

In several respects our results agree with ||l5| and . 
We find, in the case of charge neutrality, a universal I-V- 
characteristics and, over a certain range of temperature, 
a logarithmic correction to the conductance. In some 
points, however, our results are remarkably different from 
[ p5| , 16 . For the closed system our improved treatment 
of the "zero mode" allows a more precise calculation of 
the anrplitude of the In T behaviour in the conductance 
than 131. The major difference however concerns the 
low temperature saturation of dl/dV. We find that the 
scale for this saturation is set by the energy of the lowest 
diffusive mode in the system. This scale seems to be 
absent in |p^ , and the origin of this discrepancy is not 
clear to us. 

Finally, as far as the experiments are concerned, a 
InT behavior in the linear resistivity of a short metallic 
bridge together with an /-V^-characteristics which agrees 
well with the universal function (23) has been observed 
in |l^. In [|6) it has been suggested that the effect 
might be due to the Coulomb effects in a diffusive wire, 
whereas in [ [l5[ | the Coulomb correction to the tunneling 
conductance has been suggested as the explanation. Our 
work shows that also in the intermediate regime with 
both diffusive and interface resistivity the predicted I-V- 
characteristics does not change and thus agrees with the 
experimentally observed one. Furthermore we can rule 
out a purely diffusive conductor: In [ p^ the Thouless en- 
ergy, which sets the scale for the lowest diffusive mode in 
an open system and therefore the low temperature sat- 
uration of the conductance, has been estimated to be of 



the order of several Kelvin, whereas the InT is observed 
down to lOOmK. In the case with resistive interfaces on 
the other hand the energy of the lowest diffusive mode is 
reduced, 70 ~ TiD / L"^ [R^i^^ I R) < hD/L'^. From this 
consideration we conclude that in the experiment the 
diffusive resistance is considerably smaller than the in- 
terface resistance. A further hint for the importance of 
interfaces is found from the prefactor A: For the open 
system we found A « 0.25 and A = 2GiGr/{Gr + GrY 
in the tunnel limit. The experimental values |15| are be- 
tween A « 0.43 . . . 0.7, i.e. closer to the the tunnel limit 
than to the open system. In order to check these ideas 
it would be of interest to modify experimentally the re- 
sistance of the interface relative to the short bridge and 
observe both a change of the prefactor of the In T behav- 
ior and of the saturation temperature. 

In conclusion, we calculated the Coulomb interaction 
contribution to the current through structures which 
are composed of diffusive pieces and resistive interfaces. 
Our general formalism agrees with earlier studies on 
the Coulomb correction to the tunneling conductance 
I, ^ |l|l and on the Coulomb correction in diffusive 
conductors ||l|, ^ . In contrast to those earlier studies our 
formalism treats both effects on equal footing and is valid 
even far from thermal equilibrium. We concentrated on 
the zero dimensional limit valid for temperatures below 
the Thouless energy and shown that our theoretical re- 
sults provide an explanation of the experimental findings 

of El. 
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